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We study the g-dependent susceptibility x(q) of a series of quasiperiodic Ising 
models on the square lattice. Several different kinds of aperiodic sequences 
of couplings are studied, including the Fibonacci and silver-mean sequences. 
Some identities and theorems are generalized and simpler derivations are 
presented. We find that the g-dependent susceptibilities are periodic, with 
the commensurate peaks of x(q) located at the same positions as for the 
regular Ising models. Hence, incommensurate everywhere-dense peaks can 
only occur in cases with mixed ferromagnetic-antiferromagnetic interactions 
or if the underlying lattice is aperiodic. For mixed-interaction models the 
positions of the peaks depend strongly on the aperiodic sequence chosen. 

KEY WORDS: Ising model; Z-invariance; quasiperiodicity; golden ratio; 
silver mean; correlation functions; wavevector-dependent susceptibility. 

1. Introduction 

In our most recent paper,® we have studied the g-dependent susceptibihty 
x(q) for a Z-invariant ferromagnetic Ising model on Penrose tiles. (The 
is in many ways equivalent to the structure function determining diffraction 
patterns.) We have found that x(ct) is aperiodic and has incommensurate 
peaks which are everywhere dense, though only a limited number of them 
are visible at temperatures far away from the critical temperature. This is 
very different from the behavior of Fibonacci Ising models defined 
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on regular lattices, ''^'^ where x{^) is periodic and has only commensurate 
peaks located at the same positions as for the regular Ising models when the 
couplings between the spins are ferromagnetic. 

The periodicity of when the lattice is regular, is due to the fact that 
we may write 

kBTxiq.,qy) = e'('?^'+'?'''")C(/, m) (1.1) 

l,m 

where the average of the connected correlation function for two spins with 
fixed separations {l,m) is 

C{l,m) = lim [{<Jl',m'Cri'+l,m'+m) - {cri',m'){iyi'+l,m'+m)]- (1-2) 

l',m' 

in which C denotes the number of rows and columns in the lattice, so that 
£^ is the total number of spins. Since / and m are integers, it is easily seen 
from (11. ip that the g-dependent susceptibilities for such cases are periodic 
with periods 2n in and Qy. When the lattice structure is quasi-periodic, 
as in the case of the Penrose tiles studied in our previous paper, it is not 
possible to split the summation in the susceptibility in this way and x(q) is 
no longer periodic. 

In this paper, we want to examine the g-dependent susceptibility of some 
other aperiodic ferromagnetic Ising models defined on regular lattices, to 
find out if the Fibonacci Ising models are different from other more general 
aperiodic models. 

To be more specific, we consider the Z-invariant inhomogeneous Ising 
mo del ^'^'^'^'^ defined on a rectangular lattice as shown in Fig. 1, and let 
either one of the sequences of rapidities, {un)nez or {vm)mez or both, be cer- 
tain aperiodic sequences. In doing so, we shall derive a number of properties 
for these sequences which are part of the main results of this paper. 

As before, the edge interactions are parametrized by (see Fig. 2) 

sinh (2K{ui, vj)) = k sc{ui — vj, k') = cs(A + Vj — Ui, A;') , 
sinh {2K{ui,Vj)) = cs{ui — Vj, k') = ksc{\ + Vj — Ui, k'), 

(1.3) 

where A = K{k') is the complete elliptic integral of the first kind, k and 
k' = ^/l — k^ are the elliptic moduli, which are temperature variables, and 
they are the same for all couplings. 



Q-Dependent Susceptibilities of Quasiperiodic Ising Models 



3 



\ A A A 



AAA 























• 




• 

(m,n) 




{rn,n) 


• 


». 

( 




• 




• 




• 










• 

HMD 




• 




• 


^ 

( 




• 





• 




• 










• 




• 




• 























■ ^m—n+l 



^j-k-l 



■Uj-k+2 



Fig. 1. Fig. 1. The lattice of a two-dimensional Z-invariant Ising model: The rapidity 
lines on the medial graph are represented by oriented dashed lines. The positions of the 
spins are indicated by small black circles, the positions of two of the dual spins by white 
circles. 
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Fig. 2. Fig. 2. (a) Horizontal coupling Kij = K{ui,Vj); (b) Vertical coupling Kij = 

K{ui,Vj). 
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2. Quasi-Periodic Sequences 

Quasi-periodic sequences were first used — within the related context of the 
study of the specific heat of layered Ising models — by Tracy. ''^'^ Even though 
the particular sequences used by Tracy® may all be interesting, for some 
technical reasons we shall consider here only the aperiodic sequences which 
were studied by de Bruijn.® Let 

a, = i[(j + l) + v/(j + l)^ + 4], forj = 0,l,2,..., (2.4) 

such that ao = (1 + V^)/2 is the golden ratio and ai = 1 + \f2 is the silver 
mean. Define for each j a sequence {pj{n))n^i, 

Pj{n) = [7 + (n + l)/a,J - [7 + n/«,J, (2.5) 

where \_x\ is the largest integer < x, and 7 is a real number. In this paper, 7 is 
chosen such that '-f+m/aj does not equal an integer for any m. Consequently, 
the sequence in (12.51) is not changed when the floor ( [x\ ) in (12.51) is replaced 
by ceiling or roof [Ix]: smallest integer > x). For the silver mean sequence 
{pi{n))nez, we choose 7 7^ m -|- /a/2 for all integers m and /. More generally, 
it is sufficient to require that 7 is not a solution of a quadratic equation with 
integer coefficients. 

It is shown by de Bruijn® that the {pj{n))n(^z are sequences of O's and 
I's, which may also be easily shown by rewriting (12.51) as 

Pj{n) = [xn + Xn = {1 + n/aj}, (2.6) 

after decomposing j + n/aj into its integer and fractional parts, i.e., 

7 + (n + l)/aj = [7 + n/aj\ + {7 + n/aj} + l/aj, (2.7) 

{x}=x-[x\, < {x} < 1. (2.8) 

Since aj > 1 and therefore < a;„ -|- l/ctj < 2, it follows that Pj{n) = 
if < Xn + l/dj < 1 and Pj{n) = 1 if 1 < a;„ -|- l/ctj < 2. As j increases 
(and so does aj), the corresponding sequences {pj{nj)n£z contain increasing 
numbers of zeros. In fact, for fixed j, the p^-'s can be separated into blocks of 
j + 1 digits (a one followed by j zeros) or j + 2 digits (a one followed by j -|- 1 
zeros). Furthermore, it is also shown by de Bruijn® that the production 
rule of replacing each 1 in a sequence pj's by a 1 followed by j -|- 1 zeros and 
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replacing each by a 1 followed by j zeros produces a new sequence of p'j's 
of the form f l2.5p with 7—^7' and 7' = — {7}/ttj. 

All these sequences are known to be aperiodic. Thus, if we let Um = ma 
for pj{m) = 1, and Um = mb for Pj{m) = 0, then the sequence of line variables 
or rapidities (Mm)mez is related to the sequence {pj{rn))m<=:Z, and therefore is 
also quasiperiodic. For j = 0, the po{mys and the corresponding tijYi s are 
Fibonacci sequences, and this case we have studied earlier. ''^'^'^ Likewise, we 
may also associate a sequence of rapidities (fn)nez to the sequence {pj{n))n£z- 
In this way, we can construct several quasiperiodic Z-invariant Ising models 
on the square lattice. 

In order to calculate the average of the connected correlation functions, 
C{1, m) given by (11. 2p . we need to generalize a result of Tracy''^ for Fibonacci 
sequences. Tracy' ^' mentions also some other quasi-periodic sequences, for 
which — as far as we know — the corresponding theorems are not yet available. 
But we can generalize his result to general j > while simplifying his proof 
at the same time. 

Averages: Following Tracy,® we let N{n,m) be the number of I's in the 
subsequence pj{m), ■ ■ ■ ,pj{m + n — 1) which is also the number of ma's in 
the subsequence Um, - ' ' y "^m+n-i- 

Because the only allowed values of Pj{n) are either 1 or 0, the number of 
I's among these n consecutive terms of pj's is 

n-l 

N{n, m) = pj{m + i) = [7 + (m + n)/ aj\ — [7 + m/ aj\ 



i=0 



[xm + n/aj\ = [xm + {n/aj}\ + [n/aj\, (2.9) 



where Xm is defined in (12. 6p and < < 1- Since < Xm + {n/aj} < 2, 
we find 

AT/ \ i \n/aA for Xm + in/ai}<l, -.^x 

N(n,m) = <^ L / jj r / 1 ^ . (2.10) 

^ ' ^ \ [n/aj\ + 1 for Xm + {n/aj} > 1. ^ ^ 

Noting that l/aj is irrational, we find from Kronecker's theorem OS) that as 
m varies from —00 to cxo, the x^'s in (12.61) are distributed everywhere dense 
and uniformly between and 1. Thus the probability of finding an Xm with 
< 1 — {n/aj} is 1 — {n/aj}, whereas the probability of > 1 — {n/aj} 
is {n/aj}. Consequently, we have proved the following theorem: 
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Theorem 1: An infinite quasiperiodic sequence (um)mez defined by Um = ua 
if Pj(m) = 1 and Um = ub ii Pj{m) = with the p/s given by (12.51) contains 
blocks of a single ua followed by either j or j ' + 1 u^s. The number of 
Ma's among n consecutive u's is either \n/aj\ with probability 1 — {n/aj} 
or \n/aj\ + 1 with probability {n/aj}. 

We have thus generalized the result of Tracy''^ for Fibonacci sequences (with 
j = 0) to other cases (j > 0), while also simplifying the proof. 

Sequences of Three Objects: Since each pj sequence is quasiperiodic, if 
we shift a pj sequence by a certain number of digits and subtract the shifted 
sequence from the original one, the resulting sequence is also quasi-periodic, 
having three different values: 1, and —1. Moreover, as a.pj sequence consists 
of blocks of J ' + 1 digits with a one followed by j zeros or blocks of j + 2 digits 
with a one followed by j + 1 zeros, we find for j 7^ or 7^ (1 + a/5)/2, that 
two consecutive terms in the original pj sequence cannot be simultaneously 
1. Consequently, if we let 



the average number of I's (or — I's) among n consecutive numbers can be 
easily evaluated. In this paper, we restrict ourselves to the sequences (12. lip 
with j > 1. Therefore, we work out the needed probabilities next. 

As it is indeed impossible to have both pj{i + 1) and Pj{i) equal to 1, we 
find qj{i) = 1 if Pj(£ + 1) = 1; qj{i) = —1 if Pj(£) = 1, and otherwise when 
Pj{i+1) = Pj{tj = 0. Now we let Um = ua if Qjim) = 1, Um = mb if Qjifn) = 0, 
and Um = uq if Qji^) = — 1- Consequently, the sequence of rapidities Um is 
related to the sequence qj{m), and is therefore also quasiperiodic. 

Let the number NA{n,m), {N-Q{n,m) or Nc{n,m)) denote the number of 
I's, (O's or —I's) in the subsequence qj{m),- ■ ■ ,qj{m + n — l), which is also 
the number of ua, (mb or uq) in the subsequence Um+n~i- Since the 

number of I's in qj{m), ■ ■ ■ , gj(m+n— 1) is equivalent to the number of I's in 
Pj{m + 1), ■ ■ ■ ,Pj{m + n), we find 



q,{e)=pj{i + l)-p,{i), iez, 



(2.11) 



NA^n, m) 



[7 + (n+m + l)/ajj - [7 + (m + l)/ajj 
[xm+i + n/ aj\ 

f [n/aj\ for x^+i < 1 - {n/aj], 
\ \n/oij\ + 1 for Xm+i > 1 — {n/aj}. 




(2.12) 
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Fig. 3. Fig. 3. The regions of Xm where A'a(«, rn) — ln/aj\ +^ and Nc{n, m) — [«/ckjJ +1^ 
are shown for the silver-mean case j — 1. The segments where n or v ^ are indicated by 
thick white strips, while the segments where fi or v = 1 are indicated by narrow shaded 
strips. The /^-strips are below and the z/-strips on top. 



cf. (12.91) and (12.101) . Likewise, the number of — I's in the new subsequence 
qj{m), ■ ■ ■ , gj(m + n — 1) is equivalent to the number of I's in the original 
Pj{m), ■ ■ ■ ,pj{m + n — 1), and we find 

iVc(n, m) = [7 + (n + m)/aj\ - [7 + m/ajl 

= [Xm + n/ aj\ 



[n/aj\ for < 1 — {n/aj}, 
[n/aj\ + 1 for > 1 — {n/aj}. 



(2.13) 



Since the total must be n, we have 

A^b(^, m) = n — N^in, m) — Ncin, m). (2-14) 
Using (12.61) . we find 

r ,1/1 / x,^ + l/oj for + l/aj < 1, , . 

+ I ' / Ji I Xm + l/ftj-l for Xra + l/cyj > I- 

In view of the above, let us write 

Nx{n,m) = \ n/aj\ + ^, with /i = or 1, 

Nc{n,m) = [n/aj\ + u, with z/ = or 1. (2-16) 

Eqs. (12.121) . (I2.13p . and (12.151) determine the proper choices of fi and u as 
functions of Xm and {n/aj}. This is illustrated in Fig. 3 for the case j = 1; 
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the situation is qualitatively the same for all j > 1. We remind ourselves 
that the Xm defined in (12.61) is everywhere dense and uniformly distributed 
in [0,1), as m runs from —oo to oo. Consequently, the choice u = is 
seen from (I2.13P and (I2.10p to correspond to the segment in [0, 1] where the 
inequality < Xm < I — {n/aj} is satisfied, while u = 1 is given by its 
complement satisfying 1 — {n/aj} < Xm < 1- Using (12.151) . we find < 
Xjn+i < 1 — {n/aj} is equivalent to both 1 — l/a^ < < 2 — {n/aj} — 1/aj 
and < Xjn < I — {n/aj} — l/aj. Since < Xm < 1, the second inequality 
cannot be satisfied if l — {n/aj} — l/aj < 0, which is the defining condition for 
Fig. 3 (c); here /x = is the segment where l — l/aj<Xm< 2 — {n/ctj} — l/a^; 
its complement /i = 1, however, consists of two disjunct segments. Cases 
with 1 — {n/aj} — 1/aj > are shown in Fig. 3(a) and (b); here /i = 
consists of two disjunct segments satisfying 1 — l/a^ < x^ < 1 or < x.^ < 
1 — {n/aj} — 1/aj, while its complement /x = 1 is now just one segment given 
by 1 — {n/aj} — 1/aj < Xm < 1 — 1/aj. 

Let P{fi', z/'), for fi', z/' = 0, 1, denote the joint probabihty for having both 
Np^{n,m) = [n/aj\ + fi' and Nc{n,m) = \n/aj\ + u' . Then P{jji',u') is the 
total length of the intersection of the segment or segments where fi = fJ^' with 
the segment where u = v' . The results are different for the three different 
regions of {n/aj}. We find 

P(0,0) = l-2{n/a,} \ 

P(l, 0) = P(0, 1) = {n/a,} ) if {n/a,} < 1/a,, (2.17) 
P(1,1) = J 

P(0,0) = 1 - {n/aj} - 1/aj 

P(l, 0) = P(0, 1) = 1/aj ) if 1/aj < {n/aj} < 1 - 1/aj, (2.18) 
P(l, 1) = {n/a,} - 1/a, 

P(0,0) = 

P(l, 0) = P(0, 1) = 1 - {n/a,} } if {n/aj} > 1 - 1/aj. (2.19) 
P(l, 1) = 2{n/a,} - 1 

Remark: Both Theorem 1 for the two-object case and Eqs. (I2.16P through 
(12.190 for the three-object case have a refiection symmetry under the formal 
replacement n —>■ —n, (n > 0). Since aj is irrational, this means we have to 
replace {n/aj} {—n/aj} = 1 — {n/aj}, so that P(yU, z/) P(l — yU, 1 — u). 
Also, —[—n/aj\ = \n/aj] = [n/aj\ + 1, and (12.160 is to be replaced by 
Np^{—n,m) = \n/aj\ +1 — fi, Nc{—n,m) = [n/aj\ +1 — v. Therefore, we 
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indeed have that the probabihty distribution is invariant under reflections. 
It is also translationally invariant, as fl2.16l) -( l2.19l) are independent of m. 

3. Correlations 

The spin-spin correlation function in the inhomogeneous Z-invariant Ising 
model has been shown by Baxter® to depend only on the elliptic modulus k 
and the rapidity variables, m's and f 's, of rapidity lines that are sandwiched 
between the two spins. Particularly, for — Z < m < /, when the arrows of all 
these relevant rapidity lines are pointing to the same side of the line joining 
the two spins (see Fig. 1), we have — according to the rule in ref. |5] — the result 

{ci j,k(^j+l,k+m) = 92li'^j-k+l: • • • 5 Ul-m+j-k-, Vj+k-, ■ ■ ■ , Vl^m+j+k-l) , (3.20) 
{f'j,kfJ'j+l,k+m) = fi'2i(%-fc+l) • • • ) '^i-m+j-fc) • • • ; ^i+m+j+fc); (3.21) 

while for — m < I < m, when the arrows of the vertical rapidity lines and the 
arrows of the horizontal rapidity lines are pointing to opposite sides of the 
joining line, we find 

{^j,k'^j+l,k+m) = 

92mi^j-k+l-m+l, ■ ■ ■ , Uj^k, A + fj+fc, . . . , A + fj+m+j+fc-l), (3.22) 
{fJ'j,kfJ'j+l,k+m) = 

92mi'^j~k+l-m+l, ■ ■ ■ , Uj-k, A + Vj+k+1- . . . , A + Vi+m+j+k), (3.23) 

Here A = K{k') is a complete elliptic integral of the first kind. Note that the 
explicit dependence of X, g and g* on the elliptic modulus k is dropped, but 
it should still be understood to be implicitly present. Also, the fi = a* stand 
for dual spins on the dual lattice, which is at the dual temperature. 

As pointed out first by Baxter,® the universal functions g2i and (?2/ have 
"permutation symmetry" (meaning they are invariant under all permutations 
of the rapidities) and the "difference property" (which implies a translation 
invariance when shifting all the rapidities by the same amount v^^^). The 
functions (72/ ^^'^ 92i ^ > ^ "^^^ t)e obtained iteratively.''^'^'^ The final 
technical point is to explain how the averaging in (11. 2p is done. 
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3.1. Averaging 

In this paper, we shall consider quasiperiodic sequences which are either 
sequences of two objects: 

Um= I if Pjim) = 1, ( VA if Pjim) = 1, 

\ ub if Pj{m) = 0, \ fB if Pj{m) = 0, 

for fixed j > 0, or sequences of three objects: 

{Ma if qjim) = 1, 
Mb if qj{m) = 0, Vm = v, j > 1. (3.25) 
Uc if qj{m) = -1, 

To evaluate C{l,m) and C*{l,m) for \m\ < I, we use (13.201) and (13.211) . 
It is easily seen from these equations that there are / — m horizontal rapidity 
lines u and / + m vertical lines v sandwiched between the two spins. 

For the two-object sequences in (I3.24p . we find from Theorem 1 that the 
number of ma's among the l — m consecutive w's is either [sj with probability 
1 — {s} or [sJ + 1 with probability {s} where s = {l—m)/aj, while the number 
of Ma's among the / + m consecutive f 's is either [rj with probability 1 — {r} 
or [rJ + 1 with probability {r} in which r = (/ + m)/aj. Consequently, the 
averaged connected correlation function in (11.21) for |m| < / becomes 

C(/,m) = {l-{s}){l-{r})g[[s\,l-m-[s\,[r\,l + m-[r\] 

+ {1- {s}){r}g[[s\,l - m- [s\,[r\ + 1, / + m - [rJ - 1] 

+ {s}(l - {r}) g[[s\ +l,l-m- [s\ -l,[r\,l + m- [rJ] 

+ {s}{r} g[[s\ + l,l-rn - [s\ - 1, [rJ + 1, / + m - [r\ - 1] 

- {^)\ (3.26) 



where 



{l — m)/aj, r = {l + m)/aj, (3.27) 



g[m3,m2,mi,mo\ ; 

ma 



g{uA, . . . ,ua,ub, . . . ,ub,va, ■ ■ ■ , va, vb, ■ ■ ■ , vb) (3.28) 
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and (cr) = 0, as T > Tc. The averaged correlation C*{l,m) of the disorder 
variables, which is also the correlation function for T <Tc, can be obtained 
from the above equations simply by replacing g by g* and (a) by (1 — 

For the three-object sequences in fl3.25p . the numbers of ma's and uqs 
among the / — m consecutive u's are given by (12.121) and (12.131) . and the 
averaged connected correlation in (11.21) can be evaluated using (12.171) through 
dSHD- We find, for {s} < l/aj, 

C{l,m) = (l-2{s})^[[sj,/-m-2[sj, [sj,/ + m] 

+ {s}^[[sj,/-m-2[sj -1, [sj +l,/ + m] 

+ {s} g[[s\ + 1,1 — m — 2[s\ — 1, [s\,l + m] 

- (3.29) 



where 



g[m3,m2,mi,mo] 



g{uA, . . . ,ua,ub, . . . ,ub,uc, ■ . ■ ,uc,v, . . . ,v). (3.30) 
For l/aj < {s} < 1 — l/aj, we find 

C{l,m) = {1 — {s} — l/aj) g[ls\,l — m — 2[s\, [s\,l -\- m] 
+ {l/aj)g[[s\,l-m-2[s\ -1, [s\ +l,l + m] 
+ {l/aj)g[[s\ +l,/-m-2[sj - 1, [s\,l + m] 
+ {{s} -l/aj)g[[s\ +l,/-m-2[sj -2, [sj +l,/ + m] 

- (^)^ (3.31) 

whereas, for {s} > 1 — l/aj, 

C{l,m) = (l-{s})^[[sj,/-m-2[sj -l,[sj +l,/ + m] 
+ (1 - {s}) g[[s\ +1,1 -m- 2[s\ - 1, [s},l + m] 
+ {2{s} - l)g[[s\ + l,/-m-2[sj -2, [s\ +l,l + m] 

- {af. (3.32) 
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Again, the formulae for C*{l,m) are similar, cf. the discussion below (13.281) . 
Also, it is easily verified that we have the general inversion symmetry 

C{-1, -m) = C{1, m), C*{-1, -m) = C*{1, m), (3.33) 

valid for all values of / and m. Hence, we have now the results for |m| < 

To evaluate C{—m,l) and C*{—m,l) for \m\ < I, we let / —>■ —m and 
m ^ / in (13.221) and (13.231) . and find that there are I + m horizontal lines u 
and / — m vertical lines v sandwiched between the two spins. 

If Um and Vn are given by (13.241) . Theorem 1 can again be used to find 
the average number of ma's among the / + m consecutive m's and the average 
number of f a + A's among the / — m consecutive f 's. As a consequence the 
averaged connected correlation function in (11.21) for |m| < / becomes 

C{-m,l) = {l-{r})il-{s})g'[[r\,l + m-[r\,[s\,l-m-[s\] 

+ {r}il-{s})g'[[r\ + IJ + m - [r\ - 1, [s\,l - m - [s\] 

+ {1- {r}){s}g'[[r\J + m- [rj, [sj +l,/-m- [sj - 1] 

+ {r}{s}g'[[r\ + 1,1 + m - [r\ - 1, [s\ + 1, 1 - m - [s\ - 1] 

- (^)^ (3.34) 

where 

g'[m-i,m2,mi,mo\ = 

ma 1712 nil mo 



g[uK, . . . , Ma, Mb, • • • , Mb, A + Ma, . . . , A + f A, A + fB, . . . , A + Mb). 

(3.35) 

For the three-object sequences in (13.251) . the average numbers of ma's and 
Mc's among the / + m consecutive m's are given by (12.171) through (I2.19p . For 
|m| < I and {r} < l/a^, we find 

C(-m,/) = (l-2{r})^'[[rj,/ + m-2[rj, [rj,/-m] 

+ {r}^'[[rj,/ + m-2[rj - 1, [rJ + l,/-m] 

+ {r} (7[[rJ + 1, / + m — 2 [rJ — 1, [rJ , / — m] 

- W', (3.36) 
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while, for l/oj < {r} < 1 — l/aj, 

C{-m,l) = (1- {r}-l/aj)g'[[r\,l + m-2[r\,[r\,l-m] 
+ {l/aj)g'[[r\,l + m-2[s\ - 1, [rj +l,l-m] 
+ {l/aj)g'[[r\ + 1, / + m - 2 [rJ - 1, [rJ , / - m] 
+ i{r}-l/aj)g'[[r\ + 1,1 + m - 2[r\ - 2, [r\ + 1,1 - m] 
- (^)^ (3.37) 
whereas, for {r} > 1 — l/aj, 

C{-m,l) = {l-{r})g'[[r\,l + m-2[r\-l,[r\+l,l-m] 
+ (1 - {r}) g'[[r\ + 1,1 + m - 2[r\ -l,[r\,l -m] 
+ {2{r}-l)g[r\ +l,/ + m-2[rj -2, [rJ +l,l-m] 
- (t^)', (3.38) 



where 



g [1713,1712, mi, mo\ 



m3 m2 mi mo 



5f(uA, . . . , ua, ttB, . . . , ^iB, Mc, . . . , Mc, A + f , . . . , A + (3.39) 

In view of (13.331) and the discussion below (I3.28p . we have now obtained 
a complete set of formulae for C{l,m) and C*{l,m). We can use difference 
equations to obtain, by iteration, all needed g^s and (?*'s. The details of such 
calculations are in our previous work,®EO gj^f^ -yyju ^ot be presented here. 
Since the various ^ffms, m2, mi, mo]'s are obtained iteratively from g^s and 
g*'s with smaller m^'s, it is necessary to evaluate the (^[ms, m2, mi, mo]'s for 
almost all rrii such that < < even though for each fixed aj, only a 
fraction of these g^s are needed. In spite of powerful modern computers, these 
calculations are still quite time consuming. Therefore, it is more economical 
to obtain the correlations for all different j's studied in one shot. 

In the next section we shall be a little more specific. If the four above 
rapidities (ma, ub, f a, ^^b) or (ma, ub, uq, v) are chosen to be multiples of A/4 
in a certain way (details will be given later), we can use the permutation 
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property to arrange the rapidities in g and g* in descending order, and then 
use the difference property to make the smallest rapidity identically equal to 
zero. Then, all functions g and g* in ([32HD, (^M>, flsAl) and (1339]) can be 
brought to the form 

m-j 1712 "T-i rriQ 

g[m3, m2, mi, mo] = ^(f A, . . . , |A, ^A, . . . , |A, ^A, . . . , ^A, 0, . . . , 0), (3.40) 
possibly permuting the mj's. 

4. Wavevector- Dependent Susceptibility 

Since the correlation functions decay exponentially (T ^ Tc), we need to put 
all the terms that have approximately the same order of magnitude together. 
More specifically we write 

oo 

X{q., qy) = kBTxiq., qy) = C(0, 0) + 2 Si, (with C(0, 0) = 1), 

(=1 

I 

Si = [C{l,m) cos{qJ + qyTJi) + C{—m,l) cos{—qxm + Qyl)] , (4.4:1) 

m=l—l 

where 5*/ contains the correlations of the top and right edges of the square 
whose four corners are (±/, ±/). The above cosines result from the use of the 
inversion symmetry (13.331) in order to include the contributions of the other 
two edges. For T away from T^., only a few Si for / small are numerically 
significant. As T ^ T^, more and more terms need be included. This way 
the g-dependent susceptibility can now be evaluated for different cases. 

4.1. Sequences of Two Objects, Example I 

We shall first consider some quasiperiodic sequences of two objects. Let the 
sequence of rapidity lines be defined by fl3.24p with the particular values 

Ma = 3A/4, ub = 2A/4, va = A/4, vb = 0. (4.42) 

Comparing (13:281) and (13:351) with (13:401) . we find 

g[ni3,m2,mi,mo] = g[m3,m2,mi,mo], 

g'[mi, mo, mg, ms] = g[m3, ms, mi, mo], (4.43) 



Q-Dependent Susceptibilities of Quasiperiodic Ising Models 



15 



where the permutation property and the difference property® are used for 
the second identity. Now, comparing (13.261) with (13.341) . we find 

C(-m, /) = C{1, m), C*{-m, I) = C*{1, m). (4.44) 

From (14.441) . we find that the g-dependent susceptibihty must have fourfold 
rotational symmetry. 

We can calculate the g-dependent susceptibility for fixed T ^ Tc, (k ^ 1), 
to arbitrary precision using an algorithm of polynomial complexity. For 
this purpose, we use quadratic difference equations^'^'^'^'^ to numerically 
evaluate the averaged correlation functions given by (I3.26p . (I4.43P and (I4.44p 
for T > Tc, and replace g and {a) = by g* and (cr) = (1 — A;"^)^/® for 
T < Tc. We have used Maple software for this, as higher and higher precision 
arithmetic is needed closer and closer to Tc. Substituting the results into 
(14.4 ip . we obtain the g-dependent susceptibility at different temperatures. 
We shall present our results mostly in density plots to get an overview of the 
full (ga;, gy)-dependence. Our results, however, are far more accurate than 
these plots suggest. 

In Fig. 4, we show four density plots of l/x(q) for j = 0, 1, 2 or 3 
and —271 < g^, qy < 27i, at the one temperature T > Tc for which the 
above-Tc correlation length ^ ^ sEI (In the density plots, darker means a 
relatively larger value of and x = g^, y = qy.) We find that there is 

no incommensurate behavior, for all different values of aj with j > and at 
arbitrary temperature. The peaks of ^ire at the commensurate positions 
of the ordinary Ising model, i.e. {qxiQy) = (27rm, 27rn) where m and n are 
any integers. We also find that is indeed invariant under 90° rotation. 

To look at the situation more quantitatively, making sure that there are 
indeed no incommensurate peaks, we can study x(0, g) and x(g, g). We have 
plotted xiQy q) versus g for j = 0, . . . , 4 and T < Tc in Fig. 5 (a), and also for 
T > Tc in Fig. 5 (b). As j increases, there are more B type of rapidity lines. 
This in turn means more weak bonds are present in the system. Therefore, as 
j increases, the peaks in the susceptibility decrease, as shown in these plots. 
The changes are very small, however. The plots clearly show no indication 
of incommensurate peaks. The behavior of fo^' T > Tc is not much 

different from the behavior at T < Tc, except that the peaks are sharper. 

More precisely, ^ is the row correlation length of the uniform and symmetric square-lattice 
Ising model with the same value of modulus k = (cosh^(i^~^) ± [cosh^(|'C~^) ~ 1]^^^) j 
with minus for T > T^. For T < T^, we must choose plus, while ^ is then twice the actual 
row correlation length. ''^^J 
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(a) golden ratio: j = and ao = {•\/5 + l)/2 



(b) silver mean: j = 1 and ai = \/2 + 1 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-6 -4 -2 2 4 6 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-6 -4 -2 2 4 6 



(c) j = 2 and a2 = {VT3 + 3)/2 



(d) j = 3 and as = + 2 



Fig. 4. Fig. 4. Density plots of l/xilxj^y) for cases when the sequences of rapidities 
(um) and {v,n) are quasi-periodic sequences of two objects given by (|3.24p and (|4.42[) at 
k = 0.83791870 « 8), T > T^. There is no significant j-dependence. 
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(a) T <Tc: k = 4.2309029 (b) T > Tc. k = 0.2363561688 ((, = 1) 

Fig. 5. Fig. 5. Plots of xil^l) versus x = q for the cases given by (j3.24|l and (I4.42p and 
j = 0, . . . , 4. The curves for j = have the highest value at g = 0, and the peaks decrease 
in magnitude as j increases. 

We could give more density plots for different temperatures and also 
for temperatures below and above Tc. But those plots would not be much 
different from Fig. 4. We find that as T ^ Tc, the peaks of x(ct) become 
sharper. Also, the peaks of fo^' T > are sharper than those for 
T < Tc, as the correlation length above is only half in length compared to 
the one at the dual temperature below Tc.^ But it is hard to read that off 
from a density plot. 

4.2. Sequences of Two Objects, Example II 

Instead of (14.421) . we may also choose 

Mb = 3A/4, MA = 2A/4, vb = A/4, va = 0. (4.45) 
Comparing (I3.28P and (13.351) with (13.401) again, we find 

g[m2, mg, mo, mi] = ^([ms, ma, mi, mo], 

fif'[mo,mi,m2,m3] = ^([ms, m2, mi, mo]. (4.46) 

It is easily seen that f l4.44p still holds, so that still has 4- fold rotation 
symmetry. The behaviors of x{^) ^i-re essentially the same as in the previous 
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case, except that the peaks become sharper as j increases. 

4.3. Sequences of Two Objects, Example III 

If we let 

nA = 3A/4, nB = 2A/4, i^b = A/4, t;A = 0, (4.47) 

then 

g[ms, 1712, mo, mi] = g[m3, ma, mi, mo], 

g'[mo, mi, mg, ms] = g[m3, ma, mi, mo]. (4.48) 

Consequently, (14.441) no longer holds. As a result, x{^) behaves more like that 
of the rectangular Ising lattice, which is not invariant under 90° rotations, 
but still has only commensurate peaks. Density plots are shown in Fig. 6. 

4.4. Sequences of Three Objects, Example IV 

We now let the sequence of rapidity lines be defined by (I3.25P and let 

Ma = 3A/4, Mb = 2A/4, uq = A/4, v = 0. (4.49) 
Comparing (13.301) and (13.391) with (13.401) . we obtain 

^[ma, ma, mi, mo] = ^([mg, ma, mi, mo], 

^'[mo, mg, ma, mi] = g[m3, ma, mi, mo]. (4.50) 

We evaluate the x(q) i^i (I4.4ip by substituting this equation into (13.291) . 
(Km . ([S32D, (^M, flX^ and (Km . The probabilities for the three-object 
sequence given by (12.171) through (12.191) are quite complicated. Nevertheless, 
we find similar behavior for all different j's and temperatures. There is 
no incommensurate behavior — the peak of the susceptibility x(q) is at the 
commensurate position of the ordinary Ising model, {qxiQy) = (0,0), and 
repeated periodically with periods 27r. 

In Fig. 7, four density plots are presented for T < Tc at k = 1.1934332 
and for j = 1, . . . , 4. We again find that x(q) decreases as j increases. Since 
only the («„) sequence is aperiodic, the distortion due to the quasiperiodicity 
on x(ct) is along the diagonal. In this particular case, we find that the two 
diagonals are the symmetry axes of x(q)- 
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(a) golden ratio: j = and ao = (\/5 + l)/2 
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(b) silver mean: j = 1 and ai = \/2 + 1 



f • 
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4- 
2- 
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-2- 
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-6 -4 -2 2 4 6 



(c) j = 2 and 02 = (VIS + 3)/2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1 

6-4-20 2 4 6 



(d) j = 3 and 03 = ^5 + 2 



Fig. 6. Fig. 6. Density plots of l/x{Qx,qy) for the cases defined by (|3.24p and (|4.47p at 
T >Tc, k = 0.49127583 « 2). The susceptibility is like the one of the rectangular Ising 
model, as can be seen with some effort, with peaks still at the commensurate positions. 
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Fig. 7. Fig. 7. Density plots of ^/xilxT^y) for cases when only {um) is a quasi-periodic 
sequence. They are given by ((X^ and at T < Tc, fc = 1.1934332. Again the peaks 

are only at the commensurate positions. They are elongated in a diagonal direction. 
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4.5. Sequences of Three Objects, Example V 



If instead of fICTD . we let 



MA = 3A/4, uc = 2A/4, ub = A/4, v = 0, 



(4.51) 



then 



g[m3, mi, 1712, mo\ = g[ms, m2, mi, mo\. 



g'[mo,m3,mi,m2] = g[m3,m2,mi,mo]. 



(4.52) 



The resulting g-dependent susceptibility is less symmetric. Four density plots 
at T < Tc, k = 1.1934332, are shown in Fig. 8 for j = 1, ... ,4. Again we 
only find commensurate peaks. 

5. A Mixed Case 

We have examined quasiperiodic Ising lattices on a square lattice, whose 
interactions are quasiperiodic and ferromagnetic, and we have found very 
similar commensurate behaviors. 

Things change dramatically if we consider mixed cases with both ferro- 
and antiferromagnetic interactions, as we already know from our previous 
work that there will be many incommensurate peaks within the unit cell 
as the temperature moves close to the critical value. ''^'^ There is one new 
aspect: The results, especially the positions of the many incommensurate 
peaks, are heavily dependent on the value of j. We shall illustrate this with 
one example based on some ideas of Section 5 of Ref. [21 where several j = 
cases have been studied. 

Unlike the ferromagnetic case, we can now construct an example starting 
from the symmetric square-lattice Ising model and flipping the signs of the 
couplings by site-dependent gauge transformations. Using Theorem 1, Eqs. 
(5.17) and (5.18) of Ref. |2]are now replaced by 




(-l)L™/".J(i_2{m/a,}), 




(5.53) 
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Fig. 8. Fig. 8. Density plots of l/x{<lx,qy) for the cases given by (|3.25|) and ()4.51|) at 
T < Tc, k = 1.1934332. Still the peaks of xilxjly) are at the commensurate positions. 
The peaks are now elongated and a slight dependence on j may be observed. 
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Choosing a model aperiodic in both diagonal directions as in Section 5.6 of 
Ref. El the averaged connected correlation function now becomes 

C^^\l,m) =0(^')(/ + m)0(^)(/-m)C|5'^(/,m), (5.54) 

with CQ^\l,m) the connected pair-correlation function of the square-lattice 
Ising model. This implies that has many incommensurate peaks within 
the unit cell, and is given by 



^ ^ Xo[(lx + 2TT{l + m + l)/aj,qy + 2Tr{l-m)/aj 



(/ + l/2)2(m + l/2)27r4 

t=— oo m=— oo ^ / / \ / / 

(5.55) 

with Xo(q) wavevector-dependent susceptibility of the regular square- 
lattice Ising model. 

Density plots are given in Fig. 9 for four cases with correlation length 
^ ^ 16. Clearly, the results depend strongly on j. The case j = 3 is the most 
different as it almost looks like the periods have been halved. This can be 
explained easily since ^-k ja^ ~ n j1. 



6. Conclusions 



From the current work and our previous papers''^'^® we can draw several 
conclusions: 

• The wavevector-dependent susceptibilities x(q) of models, whose spin 
sites are on regular lattices, are always periodic. This includes cases 
when the interactions between the spins are quasi-periodic. 

• When the interactions between spins are quasiperiodic, but strictly 
ferromagnetic, x(ct) has only commensurate peaks, with behavior very 
similar to that of the regular Ising model. 

• The g-dependent susceptibilities of models on regular periodic 
lattices can have everywhere- dense incommensurate peaks in every unit 
cell, but only for cases for which the interactions between spins are 
mixed with both ferro- and antiferromagnetic couplings present.''^® 

• When the lattice is quasiperiodic — such as a Z-invariant Ising model on 
Penrose tiles — is no longer periodic but quasiperiodic and exhibits 
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(a) golden ratio: j = 0, ai = (\/5 + l)/2 



(b) silver mean: j = 1, ai = 
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(c) j = 2, i.e. ai = (VT3 + 3)/2 



(d) J = 3, i.e. as = V5 + 2 



Fig. 9. Fig. 9. Density plots of l/xiQxiQy) for the mixed case for four values of j, {j = 
0, . . . , 3), with Qx and Qy in the interval (— tt, tt) and k = 0.915398728 • • • . Now there are 
many incommensurate peaks and their positions depend strongly on j. The principal peaks 
are at {±qj,0), (0, ±qj), with qo = 27r(l - 1/ao) = 2.39996 ■ ■ ■ , Qi = 27r/ai = 2.60258 • • • , 
q2 = 27r/a2 = 1.90239- ••, qs = 27r/a3 = 1.48325- ••. This last value is close to 7r/2, 
which is reflected in figure (d). 



Q-Dependent Susceptibilities of Quasiperiodic Ising Models 



25 



everywhere-dense incommensurate peaks, even for the case of purely 
ferromagnetic couphngs. Only few of these peaks are visible within a 
given limited area of q-space when the temperature is far away from 
the critical temperature. The number of visible peaks increases as T 
approaches Tj..® 

There are many other quasiperiodic sequences. Still we have examined a 
variety of cases and believe that the above conclusions are quite generic. 

It may be interesting to consider the g-dependent susceptibility x{^) of 
the Z-invariant Ising model on the labyrinth''^'^'^ for which the distances 
between the spins are also aperiodic. To obey the symmetry, the couplings of 
pairs of spins must be related to the distances between the spins. When the 
distances are equal, the corresponding couplings must be chosen to be equal. 
Since the coupling K and ^ in a Z-invariant model are related by (11. 3p . our 
preliminary efforts in this regard have not been successful, but the model 
deserves further investigation. One thing we can predict: The g-dependent 
susceptibility fo^' the Ising model on the labyrinth, can no longer be 

periodic, and its peaks should be at incommensurate positions. 

ACKNOWLEDGMENTS 

We are most thankful to Dr. M. Widom for his interest in using exactly 
solvable models to study quasicrystals. This work has been supported by 
NSF Grant PHY 01-00041. 

REFERENCES 

1. H. Au-Yang and J. H. H. Perk, Wavevector-dependent susceptibility in Z-invariant 
pentagrid Ising model, J. Stat. Phys. ???:??-?? (2006). 

2. H. Au-Yang, B.-Q. Jin, and J. H. H. Perk, Wavevector-dependent susceptibility in 
quasiperiodic Ising models, J. Stat. Phys. 102:501-543 (2001). 

3. H. Au-Yang and J. H. H. Perk, Wavevector-dependent susceptibility in aperiodic 
planar Ising models, in MathPhys Odyssey 2001: Integrable Models and Beyond, M. 
Kashiwara and T. Miwa, eds. (Birkhauser, Boston, 2002), pp. 1-21. 

4. R. J. Baxter, Solvable eight vertex model on an arbitrary planar lattice, Phil. Trans. 
R. Soc. Land. A 289:315-346 (1978). 

5. H. Au-Yang and J. H. H. Perk, Critical correlations in a Z-invariant inhomogeneous 
Ising model, Physica A 144:44-104 (1987). 

6. J. H. H. Perk, Quadratic identities for Ising correlations, Phys. Lett. A 79:3-5 (1980). 

7. C. A. Tracy, Universality class of a Fibonacci Ising model, J. Stat. Phys. 51:481-490 
(1988). 



26 



Au-Yang and Perk 



8. C. A. Tracy, Universality classes of some aperiodic Ising models, J. Phys. A 11:L603- 
L605 (1988). 

9. N. G. de Bruijn, Sequences of zeros and ones generated by special production rules, 
Indagationes Mathematicae 84:27-37 (1981). 

10. G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 4tli 
edition, (Oxford University Press, London, 1960), Ch. XXIII Kronecker's Theorem. 

11. H. Au-Yang and J. H. H. Perk, Correlation functions and susceptibility in the Z- 
invariant Ising model, in MathPhys Odyssey 2001: Integrahle Models and Beyond, 
M. Kashiwara and T. Miwa, eds. (Birkhauser, Boston, 2002), pp. 23-48. 

12. B. M. McCoy and T. T. Wu, The Two-Dimensional Ising Model, (Harvard Univ. 
Press, Cambridge, Mass., 1973). 

13. M. Baake, U. Grimm, and R. J. Baxter, A critical Ising model on the labyrinth. 
Intern. J. Mod. Phys. B 8:3579 3600 (1994). 

14. U. Grimm, M. Baake, and H. Simon, Ising spins on the labyrinth, in Proc. of the 
5th International Conference on Quasicrystals, C. Janot and R. Mosseri, eds. (World 
Scientific, Singapore, 1995), pp. 80-83. 

15. U. Grimm and M. Baake, Aperiodic Ising models, in The Mathematics of Long-Range 
Aperiodic Order, R. V. Moody, ed. (Kluwer, Dordrecht, 1997), pp. 199-237. 



